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Abstract 

We examine in detail the detectability of a signal of diffuse high energy neu- 
trinos produced in the Milky Way by the interaction of cosmic rays (CRs) with 
the interstellar medium (ISM). We show that highly inhomogeneous galactic 
CR densities arise naturally from large scale drift effects and the inhomoge- 
neous distribution of CR sources. In particular, the CR density in the Galactic 
Center region (where the ISM density is maximal) could be an order of mag- 
nitude larger than the local CR density. Hence, the expected diffuse flux of 
neutrinos from the Galactic Center region becomes enhanced and could be 
detected at energies above ~ 10 GeV in a km 3 -size neutrino telescope within 
10 years of operation. The galactic anisotropy, which is the main signature of 
this flux, allows to discriminate this signal from other possible extraterrestrial 
diffuse fluxes of high energy neutrinos. 



1 Introduction 

A new generation of neutrino telescopes is now opening an exciting new window on 
the high energy Universe. Indeed, planned or already built large neutrino detectors 
are capable of providing us with potentially revolutionary observations in the fields 
of particle physics, cosmology, astronomy and astrophysics. 

Among different detection techniques proposed or already used, a particularly 
important role is played by under-ice and under-water optical Cerenkov detectors. 
This technique was very successfully implemented in the AMANDA detector [1] at 
the South Pole and in the Lake Baikal experiment [2]. The much larger, km 3 -size 
detector IceCube [3] is currently under construction on the Antarctic ice, while other 
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detectors (ANTARES, NEMO and NESTOR) are also being built in the Mediter- 
ranean Sea. While these under-ice/water optical Cerenkov detectors are sensitive at 
roughly TeV-PeV neutrino energies, which is essentially the range of interest in this 
work, other experiments extend the range of sensitivity up to ultra-high neutrino 
energies. Indeed, the radio Cerenkov detectors (such as RICE, ANITA, GLUE and 
SALSA) are sensitive roughly in the PeV-EeV energy range, while detectors making 
use of the Earth's atmosphere as a target volume (e.g. Fly's Eye, AGASA, HiRes 
and Auger) are sensitive at roughly EeV-ZeV energies. 

Being able to escape even from very dense astrophysical environments, and to 
travel vast distances undeflected and unattenuated, neutrinos are a precious probe 
to the high energy Universe, which is complementary to both gamma rays and cos- 
mic rays. Pioneer work on high energy neutrinos produced in various astrophysical 
processes was accomplished roughly three decades ago (see e.g. Refs. [4, 5, 6]). 
Nowadays, a primary motivation for studying neutrinos above TeV energies is the 
identification of possible individual, point-like extragalactic sources such as Active 
Galactic Nuclei (AGNs) and Gamma Ray Bursters (GRBs). It has also been pro- 
posed recently that jets accompanying core-collapse supernovae would produce de- 
tectable high energy neutrino bursts [7], which would reveal a connection between 
core-collapse supernovae and long duration GRBs. Furthermore, the existence of 
different kinds of neutrino sources located within the Galaxy has been suggested 
and studied (see Ref. [8] and references therein). However, the identification of 
point sources might be difficult to achieve, since they might be too weak to produce 
unambiguous directional signals in the detectors. Hence, also the detection of the 
corresponding unresolved diffuse fluxes (i.e. the integrated fluxes produced by all 
sources) remains an important task. In order to achieve this goal, a sound knowledge 
of all other possible diffuse neutrino signals is obviously crucial. 

At low energies, it is well known that conventional neutrinos produced by the 
decays of pions and kaons due to cosmic ray (CR) interactions in the atmosphere give 
the dominating signal [9, 10] (and this is, in fact, the only type of signal identified 
so far). At energies above few hundreds of GeV, however, this component becomes 
strongly attenuated due to the fact that the parent mesons live longer and interact 
before decaying. The picture then looks more confusing, since different high energy 
neutrino fluxes of diverse origin, but possibly similar intensity, are expected to show 
up [11]. 

Besides the extragalactic fluxes mentioned above, other expected diffuse compo- 
nents of comparable (or even larger) intensity are, for instance, the so-called prompt 
fluxes arising from the decay of short-lived charmed particles produced by CR in- 
teractions in the atmosphere (see e.g. Refs. [12, 13, 14]). The evaluation of this 
prompt component requires taking into account next-to-leading-order processes in 
the charm production cross section, which is also dependent on the behavior of 
the parton distribution functions at very low values (x < 10~ 4 ) of the fraction of 
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momentum carried by the partons, well below the range measured in present accel- 
erators. The large uncertainties coming from the small-x extrapolation procedures 
are propagated on the estimated prompt neutrino fluxes and span more than an 
order of magnitude. A characteristic signature of this (roughly isotropic) signal is 
the flavor ratio v e : : v T ~ 1 : 1 : ^, clearly different from the ratio 1:1:1 expected 
in a standard scenario for any extraterrestrial signal, due to the effects of neutrino 
flavor oscillations during propagation [15]. 

Another contribution to the high energy diffuse neutrino fluxes is that resulting 
from the interactions of the CRs present all over the Galaxy with the ambient 
gas in the interstellar medium (ISM) [5, 6, 16, 17, 18]. The very low densities in 
the ISM (pism — 1/cm 3 ) imply that essentially all mesons produced decay in this 
case without suffering any attenuation, and hence the neutrino fluxes are just the 
conventional ones from n, K and \i decays. Since these neutrino fluxes depend on the 
column density of the ISM along the particular direction considered, this component 
shows a significant anisotropy, being maximal in the direction to the Galactic Center 
(GC) (and generally enhanced near the Galactic Plane), while being minimal in the 
direction orthogonal to it. 

Although it was recognized that these diffuse fluxes of Galactic neutrinos could 
compete with the other expected contributions at high energies, a careful investiga- 
tion of this signal was missing in the literature. Indeed, inhomogeneities in the CR 
density at different locations in the Galaxy were disregarded in previous investiga- 
tions [6, 14, 16, 17, 18, 19], in which variations in the ISM density were generally 
neglected as well. On the contrary, in this work it is shown that large CR density 
inhomogeneities arise due to both large-scale drift effects and the inhomogeneous 
distribution of CR sources, producing an order-of-magnitude enhancement of the CR 
density in the GC region, where the ISM density is also maximal. Hence, this effect 
enhances the fluxes expected from the GC region for the diffuse signal of neutrinos 
produced in the Milky Way. 

A detailed study of this contribution is important on its own, since its detection 
would imply the discovery of a new high energy neutrino signal, which would be 
also a new probe of cosmic rays and the interstellar medium. But, furthermore, it 
is also crucially important as a background in order to identify a real extragalactic 
flux, and hence a new probe of the high energy Universe. As already pointed out, 
measurements of the anisotropy in galactic coordinates could be used to distinguish 
between these different possibilities. 

In the next Section, we focus on the diffusion/drift scenario that describes the 
transport of CRs in the Galaxy under turbulent conditions [20, 21, 22, 23]. While 
obtaining a sound agreement with local CR observations, we show that large inho- 
mogeneities in the galactic CR density do arise as well. In Section 3, we calculate the 
relevant neutrino fluxes and corresponding event rates, showing that the enhanced 
Galactic flux produced in the GC direction would dominate at high energies above 
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~ 10 4 GeV. Despite the small event rates, the detection of the diffuse Galactic signal 
in a km 3 -size telescope such as IceCube is shown to be feasible and realistic. Finally, 
the Conclusions are presented in Section 4. 

2 The transport of cosmic rays in the Galaxy 

According to the standard picture, galactic CRs are accelerated in supernova rem- 
nants, the expanding shock wave fronts produced by the explosion of supernovae. 
This mechanism is suitable for producing high energy CR fluxes with power-law 
spectra. The different nuclear components generated in this way may have upper 
cutoffs at the energies E cZ ~ ZE C) where E c would correspond to the cutoff of 
the proton component. In the context of the supernova paradigm for the origin of 
galactic cosmic rays, this upper cutoff may range between E c ~ 10 15 and 10 17 eV, 
depending on whether the strong shocks develop in the ISM or in the stellar wind 
of the predecessor star [24]. Higher cutoffs can be achieved, for instance, in a coop- 
erative acceleration scenario, in which the collective effect of multiple stellar winds 
gives rise to an additional acceleration process [25]. For the sake of simplicity, here 
we assume a featureless production of galactic CRs (ignoring cutoffs) and study how 
the spectrum is shaped by the CR turbulent transport. Anyhow, the neutrino fluxes 
at energies above the PeV are too faint to be detectable (see the expected event rates 
in Fig. 5), and hence the high energy end of the galactic CR spectrum should not be 
relevant for the present discussion. 

In order to characterize the transport of cosmic rays in the Galaxy, we will present 
here a simplified analytical realization of the diffusion/drift scenario [20, 21, 22, 23]. 
It should be stressed that the previous investigations on this scenario were concerned 
only with the local observations of CRs, proposing suitable explanations of the CR 
spectrum and the observed knees, as well as of the composition and anisotropy 
measurements. Here we focus our attention instead on the global inhomogeneities 
of the CR population in the Galaxy arising from the inhomogeneous distribution 
of sources and from the onset of macroscopic large-scale drifts. Indeed, the main 
goal of this work is to show the relevant role played by these global inhomogeneities 
in the production of detectable neutrino fluxes from the GC region. In order to 
check the consistency and suitability of the approach used here, the model results 
for the local CR spectrum are required to agree with the corresponding experimental 
observations. 

The region of propagation of the CR particles, the Galaxy, is often assumed for 
simplicity as a cylinder of radius R and height 2H. For definiteness, here we will 
consider a flat model of the Galaxy with R = 20 kpc and H = 2 kpc. The galactic 
disk, which contains most of the visible matter, has a height of a few hundred pc 
and is embedded in a much larger structure, the galactic halo, which is thought to 
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contain large amounts of dark matter. Associated to both the disk and the halo 
there are large-scale regular magnetic fields of a typical intensity of a few fiG, which 
are oriented to a good approximation in the azimuthal direction. In the following, 
a single extended regular magnetic field B will be considered, without making any 
further distinction between disk and halo components. 

The properties of this extended regular magnetic field are not well known, al- 
though observations seem to favor an antisymmetric structure (i.e. a field with 
opposite directions below and above the galactic plane), as would correspond to an 
AO dynamo field configuration [26]. It is important to mention, however, that here 
we are adopting a simplified picture of the regular magnetic field. Indeed, while 
Zeeman splitting [27] and polarization observations at sub-millimeter wavelengths 
[28, 29] revealed very strong toroidal magnetic fields in the central molecular zone, 
other observations also indicated the presence of nonthermal radio filaments [30] 
and threads [31] in the Galactic Center region, which appear to be directed nearly 
perpendicular to the Galactic Plane. Hence, the azimuthal symmetry assumed here 
should be regarded as a simplifying assumption that captures the global character- 
istics of the regular magnetic field in the Galaxy. The macroscopic drifts, anyhow, 
being large scale currents that dominate the diffusion of cosmic rays at high ener- 
gies, are not expected to be much dependent on the details of the magnetic field 
configurations. 

Superimposed to the regular magnetic field, a random component associated 
to the turbulent interstellar plasma is known to exist, with a maximum scale of 
turbulence of order L max ~ 100 pc. Since this random field component has a strength 
comparable to the field of the regular component, namely of the order of a /xG, the 
transport of cosmic rays in the Galaxy takes place under highly turbulent conditions. 
Existing observational data [32] show that the spectrum of inhomogeneities may be 
the same for the density of the gas and for the magnetic field, and that it is close to 
a Kolmogorov spectrum. Hence, the random magnetic field can be assumed to be 
given by a power law spectrum with Fourier components giving rise to a magnetic 
energy density dE r /dk oc £r 5/3 , for k > 2n/L max . 

A relativistic particle of charge Ze propagating across a regular field B describes 
a helical trajectory of radius sin#, where tl is the Larmor radius given by 



while the pitch angle 9 is the angle between the regular field and the particle velocity. 
In the presence of a superimposed random field component, the CRs scatter off 
the magnetic field irregularities with associated scales of order ri, changing their 
pitch angle but not their velocity. This collision mechanism, known as pitch angle 
scattering, provides an effective means of isotropization of the particle trajectories 
as long as the condition r L < L max is met, which for instance corresponds to energies 
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up to few x 10 8 GeV for protons propagating in the galactic magnetic fields. Another 
process known to contribute significantly to the diffusive transport of cosmic rays 
is the so-called field line random walk, which is the braiding and mixing of the 
magnetic field lines themselves [33, 34, 35]. 

Neglecting nuclear fragmentation and energy-loss processes such as ionization 
losses and adiabatic deceleration, the transport of high energy galactic cosmic rays 
is governed by diffusion. It should be remarked that adiabatic losses arising from 
convection in an expanding flow is an important mechanism in other contexts, e.g. 
in the escape of particles from expanding supernova remnants. However, in the 
absence of strong galactic winds, this term can be safely neglected in the transport 
equations describing high energy CRs confined in the Galaxy. Thus, CR propagation 
takes place according to a diffusion equation, 

V • J = Q, with Ji = -DijVj® , (2) 

where Q describes the distribution of sources, $ is the CR differential flux and J 
the CR macroscopic current. 

The diffusion tensor is given by 

Dij = (l>\\ - £>_l) bibj + D ± 6ij + D A e ijk b k , (3) 

where b = B /|B | indicates the direction of the regular field, Sij is the Kronecker 
delta and the fully antisymmetric Levi-Civita tensor. The diagonal components 
of the diffusion tensor are Dn (that corresponds to the diffusion in the direction 
parallel to the regular magnetic field) and D± (associated to the diffusion along the 
transverse directions), which depend on short scale turbulent fluctuations, while Da 
is the antisymmetric (Hall) diffusion coefficient related to macroscopic drift effects. 

Detailed numerical investigations on the behavior of the different diffusion co- 
efficients under highly turbulent conditions were recently performed [36, 37, 38]. 
Although the diffusion orthogonal to the regular magnetic field direction is typically 
much slower than the parallel one (unless the turbulence level is very high, since 
parallel and perpendicular motions then become similar), both D\\ and D± have the 
same dependence on energy (as long as tl < L max , which is the case of interest 
for this work). Indeed, it is found that D»,D± oc E m , where m characterizes the 
spectrum of the random magnetic field energy density, given by dE r /dk oc k m ~ 2 
[37, 38, 39]. In particular, m — 1/3 for the Kolmogorov case. Instead, the antisym- 
metric diffusion coefficient has a linear dependence on energy, D A oc E [20, 39]. As 
we will discuss further below, this stronger energy dependence of D A plays a crucial 
role in the transport of high energy galactic CRs, since the knee (E k ~ 3 x 10 6 GeV) 
and the second knee (E s k — 4 x 10 8 GeV) of the spectrum can be explained as due 
to an enhancement in the escape of CRs while going from the transverse diffusion 
dominated regime (at low energies) to the Hall diffusion (drift) dominated regime 
(at high energies). 
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Besides the turbulent diffusion/drift scenario, many other models have been pro- 
posed in order to explain the CR observations around and beyond the knee region. 
For instance, it was recently proposed that most CRs at energies above ~ 100 TeV 
could be originated in GRBs and that the knee feature could result from the trans- 
port of CRs produced by a recent nearby Galactic GRB [40]. Many other models 
for the origin of the knee are reviewed and compared in Ref. [41]. Anyhow, the 
diffusion/drift scenario is among the favored models to explain the observations of 
galactic CRs [41]. It is based on very natural assumptions, it reproduces well the 
spectral shape and the composition observations, while it additionally predicts an 
anisotropy increase correlated with the knee of the spectrum [42] , which is a partic- 
ular feature that allows to single out this scenario from all other proposals. Further 
evidence to clarify the origin of the CR knees will be available from mass group 
spectral observations and from more accurate measurements of the composition and 
the anisotropies, as expected for instance in the Kascade-GRANDE experiment [43] . 

Since we are here assuming the cylindrical symmetry of the system and the 
regular field to be azimuthal, the parallel diffusion coefficient plays no role and the 
macroscopic current is simply given by 

J = -D ± V$ + D A h x V$ . (4) 

Here, Da is considered to be a positive-definite quantity (corresponding to the dif- 
fusion of positively charged particles), while any sign changes arising from possible 
changes in the orientation of the regular field are contained in b. This agrees with 
Refs. [21, 22, 38] but differs from the notation of Ref. [20]. 

In terms of cylindrical coordinates, Eq.(2) is explicitly given by 



-—(rD — \ - — (d — \ + u — + u — 

r Or \ ± dr j dz \ ± dz J r Or z dz 



$ = Q , (5) 



where the drift velocities are 



and 



u r = -«™*> (6) 
r or 

Notice that, since — ±1, a change in the orientation of the regular field introduces 
a singular contribution to the drift velocities. 

In order to deal with analytical solutions, we will consider here a simplified 
approach in which D± is spatially homogeneous, while Da oc r [20] . Let us consider 
a regular field antisymmetric with respect to the galactic plane (i.e. = sgn(z)). 
The drift velocities are 

2D 

u r = -2D A S(z) , u z = — - sgn(z) , (8) 
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i.e. the vertical drifts are directed outwards the galactic plane, while on the plane 
there is a singular radial drift velocity directed towards the Galactic Center. 

Since the sources are located mainly in the Galactic Plane, their distribution 
is considered to be given by Q = 2h s q{r)8(z)E~ 13 , where q(r) is the radial source 
distribution, h s is a nominal vertical scale length and j3 the index of the source dif- 
ferential spectrum. As said above, we will consider that galactic CRs are accelerated 
in supernova remnants (SNRs). Hence, the source distribution will be given by the 
SNR radial profile [44] 

exp ^-p s r -^r~^j ( for r ^ 3 k P c ) ' ( 9 ) 

with a s = 1.69, f3 s = 3.33 and r = 8.5 kpc. Regarding the SNR population close to 
the GC (i.e. for r < 3 kpc), we will assume here a constant value (that matches the 
SNR population at r = 3 kpc), since a non negligible SNR density is likely to exist in 
the GC region. This assumption is indeed consistent with the available observations 
[45]. 

As the Galaxy is assumed to be flat (H <C R), one can neglect the radial gradients 
as compared to the gradients in z and the solution to the diffusion equation is given 
by 

/ h \ /l \ r R/ r 

**(£; r, .) = ET> (^) (iz£__) jf Ay g (yr) , (10) 

where w = 2HDa/tD± [20]. It should be noted that this flux, $*, is a single- 
component solution (as if all CRs were, say, protons). Further below we consider 
the total CR flux, $, which takes into account all relevant nuclear contributions. 
Also notice that, as expected from the radial drift velocity pointing towards the GC, 
the density at a given radius r receives the contribution from all sources located at 
larger radii, r > vq, while the outwardly directed vertical drifts tend to remove the 
CRs from the galactic plane. In order to make explicit the dependence of w on the 
CR energy, we can conveniently rewrite it as 

/ e \ 2/3 

where E k = 3 x 10 6 GeV is the energy of the knee and w a constant parameter 
of order unity. To reproduce the local CR observations, we will assume here the 
plausible value wq = 0.85. As shown below (see Figure 2), this choice allows to 
obtain a sound agreement with experimental observations. 

Figure 1 shows the single-component CR flux given by Eq.(10), for different 
locations (r = O,r /2,r ) on the galactic plane (z = 0), as a function of the 
adimensional, rigidity-scaled parameter e = E/ZE k . Notice that $* was multiplied 
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Figure 1: Single-component CR spectra for different locations on the galactic plane, 
as indicated, as functions of the adimensional, rigidity-scaled parameter e = E/ZE k 
(where = 3 x 10 6 GeV). The normalization is taken to be unity for the low-energy 
local solution. Notice the large density enhancement produced in the Galactic Center 
region by the onset of macroscopic drift effects. 

by a factor E^ +l ^ in order to produce flat solutions at low energies. Moreover, 
the normalization is taken to be unity for the low energy, local (r = r ) solution. 
For later convenience, this generic, single-component normalized solution with zero 
slope at low energies will be called ip{e). 

As anticipated, the spectra show a change of slope Aa = 2/3 driven by the 
Da/D± ratio. At low energies the CR transport is dominated by the transverse 
diffusion, D± oc E 1 ^ 3 , while at larger energies the onset of drift effects enhances the 
escape as Da oc E, hence producing a gradually steeper spectrum, which finally 
changes slope by Aa = 2/3 in about a decade of energy. As will be shown further 
below, after summing over all CR nuclear components in the range 1 < Z < 28 (ap- 
propriately scaling the results according to the magnetic rigidity E/Z) the envelope 
of the total spectrum nicely reproduces the observed knee (which results from the 
escape of the light component, mainly protons and He nuclei) as well as the second 
knee (which is due to the escape of the heavy component, i.e. the Fe-group nuclei). 
Another very interesting feature of the CR spectra shown in Figure 1 is their signifi- 
cant spatial inhomogeneity, especially evident at high energies around and above the 
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Figure 2: Local CR spectrum in the diffusion/drift scenario considered here, com- 
pared to experimental data sets. The main contributions to the total flux, which 
correspond to nuclei of H, He and Fe, are also indicated, as well as the assumed 
extragalactic flux, given by the dotted straight line. 



(local) knee. Indeed, in the high energy, drift dominated regime, the macroscopic 
CR currents are mainly driven by the large scale regular magnetic fields. Hence, 
large accumulations (or depletions) of CRs at different locations in the Galaxy arise 
naturally by the same mechanism that is responsible for the CR first and second 
knees observed locally. This phenomenon, being a manifestation of large-scale CR 
propagation properties, is generally also dependent on the configuration of magnetic 
fields and the distribution of sources assumed. 

In order to consider the contribution of all CR nuclear species from hydrogen 
(Z = 1) to nickel (Z = 28), we will assume that the spectra produced at the 
sources have constant (i.e. energy independent) spectral indices, which can then 
be inferred from low energy direct CR measurements. In the same vein, also the 
relative contribution of different nuclear components can be inferred from low energy 
observations. Hence, recalling the generic, locally normalized flux ip defined above, 
the total CR flux is given by 

/ E \ ~ az 

= So £/s(^J 1>(E/Z) . (12) 
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Table 1: Cosmic ray fractional abundances (at the reference energy E$ = 10 3 GeV) 
and spectral indices from hydrogen to nickel. 



In this expression, $ = 3.5 x 10~ 8 GeV _1 cm _2 s _1 sr _1 is the total CR flux at the 
reference energy E , hereafter adopted as E = 10 3 GeV, fz are the fractional CR 
abundances at the same energy, and az the low energy measured spectral indices, 
which were taken from the data compiled in Refs. [46, 47] (see Table 1). Further- 
more, we will assume an extragalactic isotropic proton component given by 

$ XG = 1.3 x 10- 28 ( 1Q / GeV ) ^ GeV-WVV" 1 , (13) 

i.e. a flux similar to that considered in Ref. [14]. Notice that an extragalactic flux, 
if isotropic, permeates the Galaxy homogeneously, unless other processes different 
from diffusion (e.g. reacceleration) do occur. 

Figure 2 shows the total local differential CR flux, along with its main compo- 
nents (protons, helium and iron), compared to data from several experiments. The 
extragalactic component is also shown separately. As can be observed clearly, a nice 
agreement between model results and experimental data is achieved even in this 
simplified analytical realization of the diffusion/drift scenario, and hence it provides 
a sound basis for the calculation of diffuse galactic neutrino fluxes. 



3 The diffuse galactic neutrino flux 

During their confinement in the Galaxy, CRs might interact with the ambient gas in 
the ISM, which is a very low density, non-relativistic plasma constituted mainly by 
atomic and molecular hydrogen. The mesons produced in these interactions decay 
before losing energy in secondary interactions, giving rise to neutrinos and photons. 
Moreover, the produced muons also decay without undergoing further energy losses 
and contribute to the neutrino signal as well. 

The neutrino flux produced along a certain direction in the Galaxy is given by 

POO _ pXmax 

<t> v {E v -b,l) = dE a™ (E) Y U (E,E V ) / da; $ N (E, x) p ISM {x) , (14) 
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Table 2: Binned radial distribution of the ISM density, Eq.(16). 



where a*™ is the inelastic proton-proton cross section, Y v is the neutrino yield (i.e. 



the mean differential neutrino spectrum produced in single proton-proton collisions), 
Pism the ISM nucleon number density, and &n(E,x) the total (galactic plus extra- 
galactic) nucleon CR flux at the galactic location denoted by x. Notice that the 
second integral involves the convolution of the CR flux with the ISM density along 
the line of sight defined by the galactic coordinates (b, I). 

The total proton-proton cross section in terms of the center of mass energy 
squared, s, is parametrized by 



CT ~ 

PP — 



35.49 + 0.307 In 2 (s/28.94 GeV 2 )] mb , (15) 



which combines accelerator and very high energy CR data [48]. The inelastic cross 
section is given by a™ el = Ka pp , where K = 0.8 is the energy independent inelas- 
ticity used in this work. 

As regards the ISM density distribution, a model consistent with the description 
given in Ref. [17] was assumed, namely 

Pism(t, z) = n(r) exp(-z/0.5 kpc) , (16) 

where the binned radial profiles are given in Table 2. It should be noticed that the 
ISM matter concentration peaks at the GC, i.e. in the same region where also the 
CR density is maximally enhanced. 

In order to calculate the neutrino yield, Y v , we simulated fixed target proton- 
proton collisions by means of the PYTHIA 6.228 program [49], a general purpose 
high energy physics event generator. This program provides an accurate represen- 
tation of event properties in a wide range of reactions, and particularly on those in- 
volving strong interactions, where multihadronic final states are produced. Lacking 
an exact description, the program makes use of a combination of analytical results 
and various QCD-based models, which are then implemented to model hard sub- 
processes, initial- and final-state parton showers, remnants and underlying events, 
fragmentation and decays, etc. The decay of long-lived particles, which are consid- 
ered to be stable in the context of collider physics applications, was enabled in our 
simulations in order to include the main galactic neutrino production channels. 

Figure 3 shows the differential spectra of diffuse galactic neutrinos for different 
directions in the Galaxy. Here and throughout, fluxes are given per flavor, but 
adding neutrinos and antineutrinos. The interactions between CR particles and 
the ISM produce roughly twice as many than v e , but at their arrival at the 
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Figure 3: Differential spectra for diffuse Galactic neutrinos arriving at the Earth 
from different directions in the Milky Way. Galactic coordinates are given in de- 
grees. Here and throughout, fluxes are given per flavor, but adding neutrinos and 
antineutrinos. Due to neutrino oscillations during propagation, the flavor ratio is 
v e : : v T — 1 : 1 : 1. For the sake of comparison, it is also shown the flux that 
would be produced from the GC direction if the CR density were homogeneous over 
the whole Galaxy (i.e. just taking into account the local CR observations). 

Earth the produced neutrinos are evenly distributed among the three flavors due 
to vacuum oscillations during propagation [15]. Since, as discussed above, both 
the ISM and the CR densities are largest in the GC region, the neutrino flux is 
maximally enhanced along the GC direction. Drift effects also enhance the galactic 
anisotropies, which are indeed the main signature to distinguish this diffuse galactic 
signal from other diffuse stationary extraterrestrial neutrino fluxes. Figure 3 also 
shows the neutrino flux that would be produced in this same direction if the CR 
density were homogeneous over the whole Galaxy (i.e. ignoring the inhomogeneous 
distribution of CR sources and all propagation effects, and just taking into account 
the local CR observations). As anticipated, a large enhancement takes place at high 
neutrino energies above ~ few x 10 4 GeV. Also notice the flux increment taking 
place at very high energies for the (b = 10°, I = 0°) direction, which is due to the 
vertical drifts that remove the CR particles from the Galactic Plane. 

In order to compare the diffuse galactic signal to other high energy neutrino 
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Figure 4: Diffuse galactic signal produced in the GC direction, compared to the 
conventional atmospheric fluxes (averaged in zenith angle), the prompt atmospheric 
flux and the Waxman-Bahcall prediction for extragalactic neutrinos produced in 
GRBs [50]. Also shown is the latest AMANDA limit on the high energy neutrino 
flux [1]. The two curves for the prompt atmospheric neutrino flux indicate the 
adopted range of small- x QCD uncertainties [19], although larger prompt fluxes are 
also possible. 

fluxes, Figure 4 shows the galactic signal produced in the GC direction, together 
with the atmospheric conventional fluxes (averaged over zenith angle), the atmo- 
spheric prompt flux and the Waxman-Bahcall prediction for extragalactic neutrinos 
produced in GRBs [50]. Also shown is the latest AMANDA limit on the high energy 
neutrino flux, obtained from their shower analysis [1] and consistent with previous 
results from the BAIKAL experiment [2]. The two curves for the prompt atmo- 
spheric neutrino flux indicate the adopted range of small- x QCD uncertainties (for 
further details, see Ref. [19] and references therein). The prompt flux, however, 
could be larger than shown here, becoming indeed the dominating high energy sig- 
nal for directions outside the GC region [19]. 

According to Figure 4, the galactic flux from the GC is the dominant non- 
conventional high energy neutrino signal. As shown before, this signal is highly 
anisotropic in galactic coordinates (see Fig. 3) and the fluxes coming from outside 
the GC region are expected to fall by more than an order of magnitude. Hence, this 
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Figure 5: Integral shower event rates as a function of visible energy E vis , corre- 
sponding to the differential fluxes shown in Figure 4. The rates are calculated for 
a km 3 -size detector taking data from the region (|6| < 10°, \l\ < 45°) for 10 years. 
The line marked "AMANDA" indicates the resulting integral spectrum assuming 
an E~ 2 power law, with no upper cutoff, normalized by the AMANDA differential 
limit [1] (which is actually given over a slightly smaller energy range). 



scenario is optimistic as regards the detection of a high energy diffuse component of 
extragalactic origin. In this respect, very useful information on the galactic diffuse 
background (for directions outside the GC region) can be drawn from the detection 
of the galactic signal in the GC direction. 

Figure 5 shows the integral event rates corresponding to the fluxes given in Fig- 
ure 4, as expected for a km 3 -size detector observing the region < 10°, \l\ < 45°) 
for 10 years of data taking. We focus here on shower (cascade) event rates produced 
by downgoing neutrinos. Notice that the angular range considered is consistent 
with the ~ 20° angular resolution expected for showers in IceCube. Being located 
at the South Pole, IceCube is capable of probing the GC region by means of down- 
going neutrinos. For TeV energies, the downgoing charged-current (CC) events 
are overcome by the much larger background of downgoing atmospheric muons. 
Instead, downgoing neutrinos can be measured by the observation of electromag- 
netic/hadronic cascades produced by neutrino interactions occurring within or near 
the detector fiducial volume [51, 52]. Furthermore, the observation of shower events 
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(without lepton tracks) is particularly sensitive to the v e CC channel [19] and allows 
an effective reduction of the conventional background relative to any high energy 
signal (being either the prompt atmospheric flux, or an extraterrestrial component). 
In fact, the visible energy that can be reconstructed from a shower event, E vis , is 
roughly the same as the incoming neutrino energy for a CC v e interaction, while 
in a neutral-current (NC) interaction E vis is smaller by a factor (y) ~ 0.3 — 0.4. 
This effect, convolved with the very steeply falling incoming spectra, and together 
with the additional suppression due to the smaller NC cross sections, yields NC 
shower fluxes which are suppressed by roughly an order of magnitude relative to 
the shower fluxes arising via CC v e interactions. The shower detection channel was 
indeed recently proposed as a novel technique to isolate the prompt atmospheric 
flux (see Ref. [19]). A further analysis of the sensitivity of the shower channel in 
IceCube, as well as its comparison with the complementary up- and down-going 
muon track channels, is accomplished in Ref. [52]. 

Despite the event rates being small, Figure 5 shows that the detection of the 
diffuse galactic signal in a km 3 -size telescope such as IceCube is feasible and realistic. 
A background subtraction of the conventional flux (which can be measured with 
large statistics at lower energies) would allow the identification of this relatively 
large galactic signal at energies above ~ 10 4 GeV. 

Let us finally comment on the implications of this scenario for the observations 
of galactic diffuse gamma-rays above TeV energies, an exciting new window on the 
high energy Universe opened by the measurements by the CANGAROO, VERITAS, 

H. E.S.S. and Milagro Collaborations. Besides the hadronic (ir° decay) channel, the 
galactic diffuse photons can be produced by other processes such as high energy 
electron bremstrahlung and inverse Compton scattering with the interstellar radi- 
ation field. Very recently, the Milagro Gamma Ray Observatory reported the first 
observations of a diffuse Galactic Plane gamma-ray signal at TeV energies [53]. 

Following a procedure analogous to that described at the beginning of this Sec- 
tion, we calculated the photon differential flux produced by cosmic ray interactions 
in the ISM. In particular, the integral flux averaged over the region (\b\ < 5°, 40° < 
/ < 100°), i.e. the so-called Milagro Inner Galaxy, is given by <f) y (E > 1 TeV) = 

I. 3 x 10~ 10 cm _2 s _1 sr _1 , which is consistent with the Milagro results [53]. Further 
work on the implications of the drift-enhanced CR densities for the gamma-ray ob- 
servations above TeV energies is under progress and will appear in a forthcoming 
paper. 

4 Conclusions 

The signal of diffuse neutrinos arising from interactions between CRs and the ISM 
in the Milky Way was reconsidered here taking into account a realistic scenario for 
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the CR transport. The combined effect arising from an inhomogeneous distribution 
of CR sources (which, according to the standard picture, are assumed as SNRs) and 
from the onset of CR drifts (which could be also responsible for the observed knees 
in the local CR spectrum) is shown to produce large inhomogeneities in the galactic 
CR density. In particular, this phenomenon could drive a large accumulation of 
CRs in the GC region, where the ISM density is also largest, thus enhancing the 
fluxes expected for the diffuse galactic signal of neutrinos. From our estimates for 
the event rates expected in a large, km 3 -size detector such as IceCube, the detection 
of this signal turns out to be quite plausible. These prospects of detection would be 
even better for larger experimental facilities, as e.g. in the proposed IceCube Plus 
project [54]. 

This neutrino flux would provide important information on the galactic CR puz- 
zle and on related aspects of great astrophysical interest: the nature and the dis- 
tribution of the CR sources, the structure and intensity of the galactic magnetic 
fields, the matter distribution of the ISM, etc. Furthermore, the detection of this 
extraterrestrial signal could be used to test for new physics beyond the standard 
model (see e.g. Refs. [55, 56, 57]). 

The distinct signature that allows to identify this diffuse galactic flux is the large 
galactic anisotropy. Indeed, the intensity from directions outside the GC region falls 
below the other expected high energy neutrino components. Hence, the detection 
of the galactic signal from the GC region could also provide very valuable infor- 
mation for understanding and estimating its impact as a background for measuring 
extragalactic neutrino fluxes in other directions. 

Furthermore, an increased population of CRs in the GC region also could have 
an impact on the production of neutrons via proton-proton interactions [58] and nu- 
clear photodisintegration due to ambient radiation fields [25, 59]. As pointed out in 
Refs. [25, 60], these neutron production and decay processes may lead to detectable 
gamma-ray and neutrino signals. Finally, let us recall that neutrons escaping from 
the GC region were also proposed as an explanation of the CR anisotropies ob- 
served by the AGASA and SUGAR experiments at EeV energies [58, 61]. However, 
these observations remain controversial, since they were not confirmed by the recent 
analysis performed with larger statistics by the Pierre Auger Collaboration [62]. 

Acknowledgments 

The author is grateful to Simona Rolli for help with the numerical calculations of 
this paper, as well as for valuable discussions. John F. Beacom, Esteban Roulet and 
Alexei Yu. Smirnov are also acknowledged for their very helpful and encouraging 
remarks. 



17 



References 

[1] Ackermann M et al, 2005 Astropart. Phys. 22 339 

[2] Aynutdinov V et al (BAIKAL Collaboration), 2003 Proc. 28th Int. Cosmic Ray 
Conf. (Tsukuba) 1353 

[3] Ahrens J et al (IceCube Collaboration), 2004 Astropart. Phys. 20 507 

[4] Berezinsky V S and Zatsepin G T, 1970 Soviet J. Nucl. Phys. 11 111 

[5] Berezinsky V S and Smirnov A Yu, 1975 Ap. Space Sci. 32 461 

[6] Stecker F W, 1979 Astrophys. J. 228 919 

[7] Ando S and Beacom J F, 2005 Phys. Rev. Lett. 95 061103 

[8] Bednarek W, Burgio G F and Montaruli T, 2005 New Astron. Rev. 49 1 

[9] Volkova L V, 1980 Yad. Fiz. 31 1510 

Volkova L V, 1980 Sov. J. Nucl. Phys. 3 1 784 (translation) 

[10] Gaisser T K and Honda M, 2002 Ann. Rev. Nucl. Part. Sci. 52 153 

[11] Learned J G and Mannheim K, 2000 Ann. Rev. Nucl. Part. Sci. 50 679 

[12] Thunman M, Ingelman G and Gondolo P, 1996 Astropart. Phys. 5 309 

[13] Martin A D, Ryskin M G and Stasto A M, 2003 Acta Phys. Polon. B 34 3273 

[14] Candia J and Roulet E, 2003 J. Cosmol. Astropart. Phys. JCAP09(2003)005 

[15] Athar H, Jezabek M and Yasuda O, 2000 Phys. Rev. D 62 103007 

[16] Domokos G, Elliott B and Kovesi-Domokos S, 1993 J. Phys. G: Nucl. Part. 
Phys. 19 899 

[17] Berezinsky V S, Gaisser T K, Halzen F and Stanev T, 1993 Astropart. Phys. 1 
281 

[18] Ingelman G and Thunman M, 1996 Preprint hep-ph/9604286 

[19] Beacom J F and Candia J, 2004 J. Cosmol. Astropart. Phys. JCAP11(2004)009 

[20] Ptuskin V S et al, 1993 Astron. Astrophys. 268 726 

[21] Candia J, Roulet E and Epele L N, 2002 J. High Energy Phys. 
JHEP12(2002)033 

18 



[22] Candia J, Mollerach S and Roulet E, 2002 J. High Energy Phys. 
JHEP12(2002)032 

[23] Candia J, Mollerach S and Roulet E, 2003 J. Cosmol. Astropart. Phys. 
JCAP05(2003)003 

[24] Biermann P L, 1995 Cosmic Winds and the Heliosphere ed J R Jokipii et al 
(Tucson, AZ: University of Arizona Press) 

[25] Anchordoqui L A et al, 2004 Phys. Lett. B 593 42 

[26] Han J L et al, 2002 Preprint astro-ph/0211197, to appear in Radio Pulsars, 
APS Conf. Ser., edited by Bailes M, Nice D and Thorsett S 

[27] Yusef-Zadeh F et al, 1996 Astrophys. J. 466 L25 

[28] Novak G et al, 2000 Astrophys. J. 529 241 

[29] Novak G et al, 2003 Astrophys. J. 583 L83 

[30] Yusef-Zadeh F, Morris M and Chance D, 1984 Nature 310 557 

[31] Lang C C, Morris M and Echevarria L, 1999 Astrophys. J. 52 6 727 

[32] Armstrong J W, Cordes J M and Rickett B J, 1981 Nature 291 561 

[33] Jokipii J R, 1966 Astrophys. J. 146 480 

[34] Jokipii J R and Parker E N, 1969 Astrophys. J. 155 777 

[35] Forman M A, Jokipii J R and Owens A J, 1974 Astrophys. J. 192 535 

[36] Giacalone J and Jokipii J R, 1999 Astrophys. J. 520 204 

[37] Casse F, Lemoine M and Pelletier G, 2002 Phys. Rev. D 65 023002 

[38] Candia J and Roulet E, 2004 J. Cosmol. Astropart. Phys. JCAP10(2004)007 

[39] Roulet E, 2004 Int. J. Mod. Phys. A 19 1133 

[40] Wick S D, Dermer C D and Atoyan A, 2004 Astropart. Phys. 21 125 

[41] Horandel J, 2004 Astropart. Phys. 21 241 

[42] T. Antoni et al (KASCADE Collaboration), 2004 Astrophys. J. 604 687 

[43] Navarra G et al, 2004 Nucl. Instrum. Meth. A 518 207 



19 



[44] Case G and Bhattacharya D, 1996 Astron. Astrophys. Suppl. Ser. 120 437 

[45] Case G and Bhattacharya D, 1998 Preprint astro-ph/9807162 

[46] Wiebel-Sooth B and Biermann P L, 1998 Astron. Astrophys. 330 389 

[47] HSrandel J R, 2003 Astropart. Phys. 19 193 

[48] Eidelman S et al, 2004 Phys. Letters B 592 1 

[49] Sjostrand T et al, 2001 Computer Phys. Commun. 135 238 

[50] Waxman E and Bahcall J, 1999 Phys. Rev. D 59 023002 

[51] Kowalski M, PhD Thesis, Humboldt-University, Berlin(2003), 
http: / / area51.berkeley.edu/manuscripts/ 

[52] Kowalski M, 2005 J. Cosmol. Astropart. Phys. JCAP05(2005)010 

[53] Atkins R et al, 2005 Preprint astro-ph/0502303 

[54] Halzen F and Hooper D, 2004 J. Cosmol. Astropart. Phys. JCAP01(2004)002 

[55] Beacom J F et al, 2003 Phys. Rev. Lett. 90, 181301 

[56] Beacom J F et al, 2004 Phys. Rev. Lett. 92, 011101 

[57] Beacom J F et al, 2004 Phys. Rev. D 69, 017303 

[58] Crocker R M et al, 2005 Astrophys. J. 622 892 

[59] Grasso D and Maccione L, 2005 Astropart. Phys. 24 273 

[60] Crocker R M, Melia F and Volkas R R, 2005 Astrophys. J. 622 L37 

[61] Bossa M, Mollerach S and Roulet E, 2003 J. Phys. G 29 1409 

[62] Letessier-Selvon A (for the Pierre Auger Collaboration), 2005 Preprint astro- 
ph/0507331, Proc. 29th Int. Cosmic Ray Conf. (Pune, 2005) 



20 



